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TRACING THE PECULIAR DARK MATTER STRUCTURE 

IN THE GALAXY CLUSTER CI 0024+17 WITH 

INTRACLUSTER STARS AND GAS 

M.J. JEE 

ABSTRACT 

Intracluster light (ICL) is believed to originate from the stars stripped from 

cluster galaxies. They are no longer gravitationally bound to individual galaxies, 

^'' but to the cluster, and their smooth distribution potentially makes them serve as 

r \ • much denser tracers of the cluster dark matter than the sparsely distributed clus- 

^ ■ ter galaxies. We present our study of the ICL in the galaxy cluster CI 0024+17 

Ph' using both Advanced Camera for Surveys (ACS) and Subaru data, where we pre- 

O ■ viously reported discovery of a ring-like dark matter structure with gravitational 

-^ ■ lensing. The ACS images provide much lower sky levels than ground-based data, 

, ^, ' and enable us to measure relative variation of surface brightness reliably. This 

analysis is repeated with the Subaru images to examine if consistent features 

are recovered despite different reduction scheme and instrumental characteris- 

tics. We find that the intracluster light profile clearly resembles the peculiar 

\0 ■ mass profile, which stops decreasing at r ~ 50" (~ 265 kpc) and slowly increases 

until it turns over at r ~ 75" (~ 397 kpc). This feature is seen in both ACS 

Q ■ and Subaru images for nearly all available passband images while the features 

O ■ are in general stronger in red filters. The consistency across different filters and 

instruments strongly rules out the possibility that the feature might come from 

any residual, uncorrected calibration errors. In addition, our re-analysis of the 

cluster X-ray data shows that the peculiar mass structure is also indicated by 

a non- negligible (3.7 a in Chandra and 2.4 a in XMM- Newton) bump in the 

intracluster gas profile when the geometric center of the dark matter ring, not 

the peak of the X-ray emission, is chosen as the center of the radial bin. The 

location of the gas ring is closer to the center by ~ 15" (~ 80 kpc), raising an 

interesting possibility that the ring-like structure is expanding and the gas ring 

is lagging behind perhaps because of the ram pressure if both features in mass 

and gas share the same dynamical origin. 
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Subject headings: gravitational lensing — dark matter — cosmology: observa- 
tions — X-rays: galaxies: clusters — galaxies: clusters: individual (CI 0024-1-17) 



1. INTRODUCTION 

Intracluster light (ICL), which was discovered by Zwicky (1951), is believed to originate 
from the stars stripped from cluster galaxies. Although details are still in dispute, a number 
of studies suggest that both the scatter of stars during the brightest cluster galaxy formation 
(e.g, Gerhard et al. 2007) and the tidal disruption of dwarf galaxies (e.g., Mihos et al. 2005) 
are among the dominant mechanisms of the ICL production. These stars are in general 
considered bound not to any individual galaxy, but to the cluster mass. If a large fraction of 
the ICL is produced during the assembly of the brightest cluster galaxies (e.g., Murante et 
al. 2004), the dynamically coUisionless property of the intracluster stars potentially allows 
us to use ICL as visible tracers of underlying dark matter at least in the central region of 
the cluster. 

Observationally, however, the quantitative study of ICL is difficult. The typical surface 
brightness of ICL is often quoted as ~ 1% or less of the night sky from the ground. At this 
faint level, instrumental systematic effects such as residual flat-fielding, scattered lights of 
bright stars, etc. become critical issues. In addition, most galaxy edges continuously blend 
into ICL, which obviously causes an ambiguity in determining where the galaxy light stops 
and ICL begins. 

In this paper, we present a study of ICL in the galaxy cluster CI 0024-1-17 at z = 0.4, 
where we recently discovered a peculiar ringlike dark matter structure (Jee et al. 2007; J07 
hereafter). The two-dimensional mass reconstruction of J07 shows that the core of the cluster 
is surrounded by a r ~ 0.4 Mpc ring-like overdense region. The feature is strongly constrained 
by coherent fiuctuation of background galaxy shapes across the r ~ 0.4 Mpc circle (weak- 
lensing) signaling the sudden change of the density slope. Of course, it is tempting to try to 
detect the peculiar dark matter structure with the cluster galaxies (e.g., Qin et al. 2008). 
However, because the density contrast of the feature with respect to the neighboring region is 
low (< 5%), the sparsely distributed galaxies cannot provide sufficient statistics to overcome 
the shot noise even if the cluster galaxy membership identification is next to perfect. This is 
the reason that here we investigate, instead, the intracluster stars, which diffusely distribute 
in the cluster potential and thus sample the underlying dark matter much more densely. In 
this investigation, the measurement accuracy is not limited by the poissonian scatter as in 
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the study using galaxies, but by the ability to control various systematics critical to precision 
surface photometry. 

We employ two sets of data: Hubble Space Telescope Advanced Camera for Surveys 
(HST/ACS) and Subaru/Suprime Cam images. Space-based imaging provides significant 
advantage over ground-based effort mainly because the sky is substantially darker (e.g., 
~ 1 and ~ 3 magnitude fainter in r and z band, respectively). Also, as the instrument is 
above the atmosphere, there is no cumbersome, time-dependent airglow effect, which often 
prevents stable calibration. However, unfortunately, there have been no ICL studies so far 
with the HST/ACS data. This is because the pipeline flatfielding accuracy ~ 1% (Mack et al. 
2002) has been considered insufficient to enable such studies. Therefore, we undertake this 
time-consuming task of producing/verifying HST/WFC flats using "blank images" before 
we proceed with the ICL measurements. One drawback in using ACS data is the small field 
of view of the instrument, which only covers a 3' x 3' region. Because the existing ACS data 
of the cluster were taken without any plan for this kind of study, the instrument pointing 
was almost fixed to the center of the cluster (with only a few tens of pixels dither to fill 
the two CCD gap), and thus a care must be taken to minimize the ICL contamination in 
determining the background level. 

Although the sky is much brighter, the Subaru/Suprime Cam images provide a few 
critical crosschecks. First, the large field of view of the Suprime Cam allows us to estimate 
the background level directly, which is important in determining the net amount of ICL and 
assessing the degree of the aforementioned ICL contamination in the HST/ACS measure- 
ments. Second, the Subaru data include the NB912 narrow band image, which Kodama et al. 
(2004) used to probe Ha emission as a measure of on-going star formation of the cluster at 
z = 0.4. The sky level in this narrow band is about 1.5 mag darker than in the broadband z' 
filter while the narrow passband gives high contrast to the Ha emission aX z = 0.4 relative 
to continuum sources at different redshifts. Third, various instrumental signatures including 
the point spread function (PSF), internal reflection, obscuration, geometric distortion, etc. 
are different from ACS. Hence, if a significant feature clearly revealed in one instrument, is 
absent in the other, this indicates that the feature might come from some residual calibration 
errors or time- dependent sources. On the contrary, detection of a consistent feature in both 
images provides evidence against these systematics. 

The hot intracluster gas also samples the underlying dark matter very densely although 
its coUisional nature biases the distribution from that of dark matter. A common isothermal 
l3 model assumes Pgas oc p^^ and tjj oc n^, where e//, and n^ are the free-free radiation 
emissivity and the electron density, respectively. Therefore, X-ray profiles in general are 
more centrally peaked than mass profiles. Also, as observed in many merging systems (Jee 
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et al. 2005; Clowe et al. 2006), some gas peaks are often significantly displaced with respect 
to dark matter peaks, which is interpreted as signaling the different physical property (e.g., 
collisional versus coUisionless nature). Therefore, a rigorous simulation is desired to fully 
address whether or not it is feasible to detect the analogous ring-like structure in the gas 
profile of CI 0024+17, given the low level of the contrast in the dark matter structure. In this 
paper, we approach the issue purely observationally by re-analyzing the existing Chandra 
and XM.M-Newton X-ray data. The current analysis is different from J07 only in that now 
we choose the geometric center of the dark matter ring as our origin in setting up radial 
bins. If there is indeed any feature in the gas profile hinting at the peculiar dark matter 
substructure, the signal is expected to be of low contrast and thus should be sensitive to the 
choice of the origin. Therefore, it is quite possible that J07 might have failed to detect the 
feature in X-ray because the center was not placed at the optimal location. 

In this paper, we adopt a Qm = 0.3, ^a = 0.3, and h = 0.7 cosmology, where the plate 
scale is 5.3 kpc/" at the redshift of the cluster z = 0.4. The quoted uncertainties are at the 
1 - CT (~ 68%) level. 



2. OBSERVATIONS AND REDUCTIONS 

2.1. ACS Data Reduction 

The cluster was observed with the Wide Field Channel (WFC) of the ACS in 2004 
November. A single pointing (~ 3'.3x ~ 3'. 3 field of view) is centered near the cluster core 
(aaooo ^ OO'' : 26"^ : 35^ (Jsooo ^ 17° : 09' : 43") with integrations of 6435 s, 5072 s, 5072 s, 
8971 s, 10144 s, and 16328 s in the F435W, F475W, F555W, F625W, F775W, and F850LP 
filters, respectively. The low level CCD processing was carried out using the STScI standard 
ACS calibration pipeline (CALACS; Hack et al. 2003) with some important modifications 
as follows. CALACS applies the LP-fiats derived from both the laboratory before launch 
(Bohlin et al. 2001) and the 47 Tuc stars on orbit (Mack et al. 2002). We find that 
these pipeline fiats deviate from sky fiats at the <0.5% level (see §2.3.11 for details). This 
level of inaccuracy should be corrected prior to ICL measurement although the sky level 
of HST images is 1-3 mags fainter. We applied these residual sky fiats to these CALACS- 
processed FLT images. Careful inspection of these results reveal that the amplifier-dependent 
(quadrant-to-quadrant variation) bias at the level of ~ 1 DN were present in some frames. 
We fix these quadrant-to-quadrant variation by comparing the medians of the two adjacent 
strips (~ 100 pixel width) evaluated after masking out astronomical objects and cosmic rays. 
In principle, these additive biases are supposed to be corrected before the application of the 
flats, which are multiplicative. However, as the bias is small and the flats are already close 
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to unity, the errors arising from reversing the order of these two operations are neghgible. 

We measured sky levels in each frame through iterative sigma-clipping and subtracted 
the values before stacking. Table 1 summarizes mean sky levels for each filter. The final 
high-level processing involving geometric distortion correction and cosmic-ray removal was 
performed using the "apsis" pipeline (Blakeslee et al. 2003). We used nearest neighbor 
interpolation in drizzling (Fruchter and Hook 2002). This is different from the procedure of 
J07, where the LanczosS (windowed sine function) kernel is employed. The LanczosS kernel 
provides a sharp point spread function, and thus suitable for weak-lensing studies. However, 
this kernel correlates neighboring pixels and artificially reduces noise fluctuation (~ 8%) 
although the degree of the correlation is found to be much lower than other commonly 
selected kernels such as square, gaussian, etc (Mei et al. 2005). Our choice of the nearest 
neighbor interpolation is to prevent any potential distortion of the intrinsic noise power 
spectrum due to drizzling. 



2.2. Subaru Data Reduction 

We retrieved the archival Subaru images of the cluster from SMOKAJ^. The data were 
taken in 2002 September 7 with the prime-focus camera (Suprime-Cam; Miyazaki et al. 
2002) in the NB912 (Ae//=9139 A, FWHM=134 A), B, R^, and z' filters with integrations of 
10,800 s, 3,600 s, 5,280 s, and 1,980 s, respectively. The number of visits per filter and the 
sky level variation are summarized in Table 2. The Suprime Cam focal plane is tiled with 
2x5 CCDs covering an area of 30' x 27'. Because we are only interested in the central region 
of the cluster, we only used the central six CCDs, which still cover a significantly large area 
(20' X 27'). 

The Suprime-Cam data reduction software SDFRED (Yagi et al. 2002; Ouchi et al. 
2004) was used to subtract overscan, mask out regions affected by Acquisition and Guide 
(AG) probe and bad pixels, perform fiatfielding (see l2.3.2l for details of our sky fiat creation), 
and correct geometric distortion. The sky subtraction procedure of SDFRED removes small- 
scale sky variations, which undesirably distorts the ICL structure and the large scale PSF 
profiles. Thus, we replace this step with our manual procedure summarized as follows. Each 
flatfielded, 20' x 27' frame shows < 1% sky gradient. This type of residual gradient is 
common and is attributed to night-by-night atmospheric effects rather than instrumental 
fiat changes. In general, the pattern is monotonic and can be removed by modeling the 
gradient with a low-order polynomial plane. In this procedure, a care must be taken to mask 
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out extended PSF wings of bright stars, galactic cirrus, and the central region of the cluster 
(3' X 3'), as well as resolved astronomical objects. Masking out PSF wings of bright stars 
was not difficult for several m^ = 10 ~ 12 mag stars inside the image. However, there is 
one rrir ^ 7.5 mag star outside (~ V away) the camera's field of view, which causes severe 
scattered light in the southwestern corner. Because of the large scale (~ 4') of the feature, 
it is difficult to disentangle this scattered diffuse hght from the sky gradient. Therefore, we 
chose to model/subtract this scattered light and the residual sky gradient at the same time by 
fitting a 4*^* order polynomial plane to the entire field. It is possible that any imperfection in 
this procedure increases a large-scale sky level fiuctuation beyond the fiatfielding accuracy. 
Therefore, we examined sky levels at random locations, and verified that the variation is 
consistent with what we expect from the fiatfielding errors. 

After we subtracted the mode from each frame, we determined the spatial alignment and 
fiux scaling between frames using bright non-saturated stars. This information was provided 
to SDFRED to create the final mosaic images. Zeropoints were evaluated by comparing 
photometry with that measured from the ACS images. We verified that these zeropoints are 
consistent with the values independently determined by Kodama et al. (2004), who kindly 
provided us with their resultqj. 

We display in the top panel of Figure [T]the color-composite image of the cluster from the 
Subaru mosaic images. The image shows the central 10' x 10' field with the blue, green, and 
red intensities representing the B, R^, and z' fiuxes. White circles denote the spectroscopic 
cluster members (0.37 < z < 0.42) that we obtained from the publicly available Moran et 
al. (2005) catalog. The yellow line delineates the 3' x 3' ACS field, which is also separately 
shown on the lower left corner. The lower right panel displaces the diffuse light of the cluster 
measured from the ACS F625W image after objects are subtracted/masked out(see [2l4|) . 



2.3. Flat-Fielding 

2.3.1. ACS Flat-fieldmg 

The ACS fiats used by the current pipeline of the STScI, called LP-fiats, were derived by 
combining the factory-measured P-fiats before launch (Bohlin et al. 2001) and the infiight 
L-fiats (Mack et al. 2002), which correct the pixel-to-pixel and the low-frequency varia- 
tions, respectively. The L-fiats were created directly for the filters F435W, F555W, F606W, 



^Kodama et al. (2004) used Vega-based magnitude system while in the current paper all the magnitudes 
are given in the AB system 
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F775W, F814W, and F850LP by observing stars in 47 Tuc, and indirectly for the remaining 
filters by interpolating the direct measurement results. Mack et al. (2005) report that the 
LP-flats are expected to reduce field-dependent sensitivity variation down to <1% for the 
first six filters, and <2-3% for the interpolated filters. Pavlovsky et al. (2006) mention that 
their investigation of the sky flats constructed from GOODS images for the F606W, F775W, 
and F850LP filters are in good (< 2%) agreement with the pipeline LP-fiats. 

The ~ 2% flatfielding error (i.e., the reported difference between the LP-fiats and the 
STScI sky fiats), if present on a large scale, is of concern for the investigation of the ICL 
features that this paper studies; for example, this will limit our ability to perform surface 
photometry to the fi < 27mag arcsec"^ regime in the F625W filter. Because the sky fiats 
of the STScI are not yet pubhcly available, it is impossible to determine the pattern of 
the deviation and thus the impact of the reported inaccuracy on our ICL measurement. 
Therefore, we decided to independently create sky flats for all six ACS filters used for the 
CI 0024+17 observations. 

We collected 90 ~ 300 blank sky images for each filter from the STScI archive. We 
avoided images, which have large extended sources (e.g., nearby galaxies), many bright stars 
(e.g., globular clusters), or diffuse light (e.g., ICL in galaxy clusters, galactic cirrus, etc.). 
The selected images were manually examined, and discarded if any noticeable sky gradient is 
present. It is well known that quadrant-to-quadrant variation (1 —3% of background counts) 
is present in WFC images (see Sirianni et al. 2002 for description of the problem). Unless 
corrected for, these features are clearly visible in the final sky fiat because the pattern is not 
completely random. We fix these quadrant-to-quadrant variation by comparing the medians 
of the two adjacent strips (~ 100 pixel width) after masking out astronomical objects and 
cosmic rays. Then, we median-smooth each image with a box size of 32 pixels, and normalize 
the image by dividing it by the mode of the image. The final, median-stacked image of 
these normalized, median-smoothed images shows the residual flatfielding error because we 
perform the task on FLT files that have already been applied the pipeline LP-fiats. We show 
in Figure |2] these residual sky flatqj. Table 3 lists the number of images that we used for the 
creation of sky flats for each filter, the deviation from the pipeline flat, and the accuracy of 
the residual sky flat that we estimate by bootstrapping. The < 0.1% accuracy in flatfielding 
implies that we can probe the surface brightness limit down to the /i ~ 30mag arcsec"^ level 
in nearly all ACS filters if flatfielding is the only dominant source of errors. 

Our sky flats confirm the claim of Mack et al. (2005) and Pavlosky et al. (2006) on the 
< 2% accuracy of the pipeline LP-flats. The most notable large scale structure is the donut- 
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like pattern particularly clear in red filters (the difference in residual between the center of 
the donut pattern and the trough (~ 1' away) is ~ 0.5% and ~ 2% in F775W and F850LP, 
respectively). This residual feature is also mentioned by Pavlovsky et al. (2006), and they 
suggest that this might be due to the difference in spectrum between 47 Tuc stars and the 
sky. While we agree that the difference in color might be a plausible source of the residual, 
we expect the residual to be still observable even if there is no color difference between 47 
Tuc stars and the sky. The sky flats measure both sensitivity and projected pixel area (i.e., 
due to geometric distortion) effects whereas the latter is difficult to measure without bias 
using stellar photometry alone unless the PSF is sufficiently oversampledjj 

On the center of the donut-like feature in F435W and F475W, there is a strong indication 
that the pipeline flats over-correct the sensitivity in this region. Comparing these residual 
sky flats with the pipeline flats also reveals that there are a few dust-moats (e.g., bottom 
of F555W and center of F625W) that were not included in the pipeline LP-flats. Finally, 
it is clear that a gridlike pattern exists in all six filters shown here. The feature is most 
(least) obvious in the F850LP (F625W) filter. The average distance between grids is ~ 60 
pixels. Because this checkerboard pattern does not change in size and location as we vary 
smoothing kernels and sizes, we conclude that the pattern is not an artifact of our 32- 
pixel median smoothing. The exact cause of the checkerboard pattern has not been known 
yet. Nevertheless, we suspect that this particular residual pattern might be the remnant of 
the ACS CCD fabrication process that is not corrected by the P-flat. The peak-to- valley 
variation of the feature amounts to ~ 1% in F850LP. 



2.3.2. Subaru Flat-fielding 

We retrieved blank images taken during the 2002 August-September period from the 
Subaru archive. As in the case for ACS, we manually examined and discarded frames that are 
not adequate for flatfield generation (e.g., presence of diffuse emission, extremely saturated 
stars, crowded stellar fields, etc.). Bad pixels and the vignetted areas by the Acquisition 
and Guide (AG) Probe on the top five CCDs were masked out in the remaining frames. 
Astronomical objects were identified and masked out by searching for 5 or greater continuous 
pixels above 1.5 times the sky rms. After median-smoothing each frame with a box size of 
16 X 16 pixels, we normalized the image using its mode, and then median-stacked all frames 



^Imagine an extreme case where the PSF is severely undersampled and the stellar profile is confined to 
a single pixel. Then, the total flux in each CCD pixel only informs us of the sensitivity, not the projected 
pixel size. 
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to create preliminary sky flats. These preliminary sky flats allow us to reflne the evaluation of 
the previous modes and the detection of astronomical objects to be masked out. Therefore, 
we obtained the flnal sky flats by iteratively evaluating the modes, detecting objects, and 
creating flats as done by Morrison et al. (1997). One important modiflcation to the Morrison 
et al. (1997) method is that we removed low-frequency residual sky gradients (< 1%) by 
fltting second-order polynomial planes. As noted by many authors (e.g.,Feldmeier et al. 
2002), each frame possesses non-negligible sky gradient due to many atmospheric effects 
even if a perfect sky flat is applied. This removal of the residual sky gradient further reduces 
the width of the sky distribution within a frame, which also helps us to mask out faint 
astronomical objects more efficiently. After a few iterations, we obtained converged master 
sky flats. 

The accuracy of these flnal sky flats were not limited by photon statistics thanks to the 
large light-collecting power of the Subaru telescope, but by the aforementioned large-scale 
sky variations in input frames. Therefore, we estimated the accuracy of the sky flat by 
bootstrap-resampling the input frames. The estimated error for the six central CCDs that 
we keep for the cluster reduction is on average ~ 0.07% (see Table 4 for the individual fllters). 
The 0.07% accuracy in flatflelding, if this is the sole source of systematics, allows us to probe 
the surface brightness ~ 8 mag deeper than the sky level. For the B fllter, this surface 
brightness limit is ~ 30.5 mag arcsec"^, similar to the value in ACS data whereas for the z' 
fllter, this limit implies that we can reach down to the fi ~ 27.5mag arcsec"^ level because 
of the high sky level (/i ~ 19.2 mag arcsec"^). However, it is important to remember that 
the quoted accuracy represents the large-scale error within the entire fleld (20' x 27'). While 
this large-scale error is a major source of uncertainty in our determination of the background 
level (thus limiting the accuracy in absolute ICL level measurement), the flatflelding error 
affecting the relative signiflcance of the ICL proflle is the uncertainty of the flats within the 
central 3' x 3' region. As this area occupies only ~ 1.7% of the total 20' x 27' area, the flat 
accuracy within the region is signiflcantly (a factor of two) better than the accuracy across 
the fleld. In addition, the dithering of the observation improves the flatness further, turning 
systematics into statistics. Because the current paper is focused on the structure of the ICL 
profile within this region, we distinguish these relative errors from the absolute errors when 
necessary. 

For a sanity check, we compare these sky flats with dome flats. For the B, z', and NB912 
fllters, high S/N dome flats were taken along with the cluster CI 0024+17 on the same night 
(2002 September 7). For the Re fllter, the closest dome flats in time were taken on 2002 
September 4. Again, we limited our analysis to the central six CCDs. We found an rms 
difference of ~ 0.2% on average, which suggests that the Suprime-Cam flats are time-stable 
at least over this two-month period. Furthermore, we note that this ~ 0.2% discrepancy is 
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dominated by a large scale gradient, and we suspect that this is caused by the non- uniformity 
of the dome screen. When this gradient is removed by second-order polynomial modeling, 
the agreement becomes < 0.1%. 



2.4. Object Detection and Masking 

As it is impossible to model the surface brightness profile of individual galaxies accu- 
rately and to subtract the contribution without introducing biases in our ICL study, we 
choose to mask out galaxies. In order to obtain consistent masking regions across filters, 
we used a single detection image for each instrument. For the ACS data set, this image is 
automatically generated by apsis from weight-averaging all six filter images. For the Subaru 
data, the i?c-band image is significantly deeper than the rest, and thus this Re image is 
chosen to detect astronomical objects. We detect objects via SExtractor by searching for at 
least 5 connected pixels above 1.5 times sky fluctuation. 

Because galaxy light continuously blends into background light, it is important to define 
the size of the masking area very carefully. An ideal masking size is to minimize the impact of 
the diffuse wing of objects while still leaving a sufficient number of background pixels usable 
for the ICL measurement. Obviously, an ellipse defined by SExtractor's semi-major and 
minor axes does not sufficiently mask out diffuse wings of galaxy light. Hence, we enlarged 
the masking size by integer multiples of these axes (we refer to these as "1 x", "2 x", "3 x", 
etc.), and examined where our ICL profile starts to converge. We observe that the axes 
of the masking ellipse should be at least three times the value given by SExtractor. This 
empirical finding is valid for both ACS and Subaru images. In Figure [3], we illustrate how 
much our determination of the location of the Gaussian peak, which we use in this paper as 
an indicator of the ICL level, depends on the masking size. In this example, the (skewed) 
Gaussian curves are the histograms of the pixel values in the F775W image at r = 76" — 84". 
With insufficient masking, the Gaussian curves are skewed because the object light spilt 
outside masking apertures contributes to the bright ends. The suitable masking size can be 
determined by examining either the skewness of the Gaussian or the shift of the centroid 
between subsequent masking schemes. In this study we adopt the latter as our guiding tool 
and find that a subsequent change in the centroid (dashed line) is at the /x ^ 30 mag arcsec"^ 
level if the axes of the masking ellipse are three times or larger than the semi-major and 
-minor axes given by SEXtractor. When we choose the segmentation map (isophotal area) 
produced by SExtractor (with the 1.5 a detection threshold) instead, the result lies between 
the 3 X and 4 x masking schemes (brown). We note that alternatively some authors set a 
very low threshold to detect astronomical objects and use the resulting segmentation map to 



- 11 - 



define the masking areas. We report tliat lowering the threshold from 1.5 a to 1.0 a would 
produce a result close to our "5 x" masking case (blue), which we conservatively choose in 
this paper. 

As for stars, we masked only the cores of the stars and subtracted the PSF wings from 
the image. We used saturated stars to sample the radial profile of the PSF wing for both ACS 
and Subaru. In the STScI archive, we found the saturated images of the star HD39060 (Beta 
Pictoris, PLPaul Kalas) for the F435W, F606W, and F814W filters. The star is located at 
the center of WFC, and this allows us to obtain the PSF profile up to ~ 150". The filter-to- 
filter variation is up to ~ 1 mag arcsec"^ at r = 100" when the profile is normalized using 
the surface brightness at small radii. We use extrapolation to model the F850LP PSF wing 
and interpolation to model the profile for the rest F475W, F555W, F625W, and F775W 
filters. For the Subaru/Suprime-Cam data, we used the brightest (m^ ~ 8.6) star in the 
CI 0024+17 image ~ T away from the cluster center. We display these ACS and Subaru 
PSF profiles in Figure HI 

Within the ACS field, we have only 8 moderately bright stars in the magnitude range 
of m^ = 16 — 19. We masked out the r = 7.5" circular area, which covers the /i/j ~ 
27 — 30 mag arcsec"^ region for these stars. The region outside this circle is PSF-subtracted, 
and we estimate that the residual error is yU/j > 31 mag arcsec"^. As can be seen in the 
top panel of Figure [H outside the ACS field (within ~ 2' from the field boundary) there are 
four very bright (m^ = 11 — 13) stars. Among these, the PSF wing of the one in the west 
(a, 5) = (00 : 26 : 24.25, 17 : 08 : 18.0) is non-negligibly affecting the southwest corner of the 
ACS image. Based on the magnitude of the star (m^ ~ 11), we expect the surface brightness 
of the PSF wing at the western edge to be fin ~ 28 mag arcsec"^. The residual error after 
PSF subtraction is estimated to be f^R > 31 mag arcsec"^, and thus should not affect our 
ICL measurement. 

The lower right panel of Figure [1] displays the diffuse light of the cluster in the ACS 
F625W image. The galaxies and the stars are masked out/subtracted as described above. 
The image shown here is obtained after median-smoothed with a box size of ~ 3" x 3". The 
color bar represents the intensity on a linear scale. The ICL distribution of CL0024+17 is 
somewhat asymmetric and extends northwest toward the secondary galaxy number density 
peak, which lies just outside the ACS field. 

For Subaru images, our PSF model extends to r ~ 5', and we subtracted all niji < 17 
mag stars while masking out the central region with magnitude-dependent radius. In the 
region overlapping the ACS field, a care was taken in choosing masking radius so that 
the estimated residual surface brightness error outside the masking aperture is below the 
~ 30 mag arcsec"^ level. 
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2.5. Background Level Determination 

The accuracy in the estimation of background sky level affects the fidelity of the ICL 
profile at large radii, and is one of the most significant sources of systematic errors in the 
current analysis. 

For Subaru data, after we removed residual sky gradients and bright stellar profiles 
( §2.21 and §2.4p . we determined the background level from the 4' < r < 10' annulus. The 
inner radius of the annulus was determined by creating the radial profile of the surface 
brightness and then locating the region where the profile starts to flatten. The uncertainty 
of these background measurements should be dominated by large scale flatflelding errors and 
residual sky gradients because the statical noise is only at the ~ 10^^ % level of the sky. We 
estimated the errors in the background level measurement by subsampling the sky within 
the 4' < r < 10' annulus. Table 5 displays these values in terms of fraction of the sky and 
the corresponding surface brightness level. 

The small fleld of view of the ACS image does not allow us to obtain the background 
level far from the cluster center. If the ICL level within the ACS field is still non-negligible, 
the background level directly measured within the ACS images biases the ICL level artifi- 
cially low. We considered performing photometric transformation from Subaru to ACS and 
estimating the corresponding surface brightness level of the ICL (thus indirectly determines 
the background levels) in ACS. However, it is difficult to prove that the two instruments in 
very different environments allow us to match the background level in this way. Therefore, 
we chose to estimate the background level still within the ACS field while minimizing the 
ICL contamination in the following way. 

Inspection of the two-dimensional images of the CI 0024+17 diffuse light reveals that 
the two-dimensional ICL distribution is highly asymmetric around the center. In case of 
F625W (e.g., lower right panel of Figure [1]) the surface brightness in the northwestern corner 
is ;U ^ 27 mag arcsec"^ whereas it is ;U > 31 mag arcsec"^ in both the southeastern and 
the southwestern corners. Because this high S/N pattern is also supported by other images 
of both ACS and Subaru, the two-dimensional feature is believed to represent the intrinsic 
distribution of the ICL. In addition, the spatial distribution of the cluster galaxies is similar 
to this two-dimensional ICL map (top panel of Figured]). Therefore, we argue that if we 
measure the background level at r > 80" from the ACS image while excluding the northwest 
region, we can minimize the influence of the ICL. This argument is further supported by our 
experiment with the Subaru image, which shows that the background level measured in this 
way (i.e., cropping the 3' x 3' region to simulate the ACS fleld) is in good agreement with 
the result measured from the 4' < r < 10' annulus within the /i ~ 30 mag arcsec"^ level. 
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2.6. X-ray Data Reduction 

For the Chandra X-ray analysis, we followed the procedures described in J07, using 
the Chandra Interactive Analysis of Observations (CIAO) software version 3.3 and the 
Calibration Database (CALDB) version 3.2. The XMM-iVewton data (taken on January 
2001 for a total integration of 52.1 ks, 52.1 ks, and 48.3 ks for MOSl, M0S2, and PN) 
were retrieved from the XMM-A^ewton Science Operations Centrqj and processed with the 
Science Analysis Software (SAS) version 7.1.2. We combined the three instrument data 
and applied the exposure map to obtain an exposure-corrected image. We detected point 
sources in the Chandra and used the results to mask both the Chandra and XMM-A^ewton 
exposure-corrected images before we measured the final X-ray surface brightness profile. 



3. ICL ANALYSIS 

3.1. Measurement of the Radial Profile 

We adopt the statistical approach of Uson et al. (1991) in measuring the radial ICL 
profile of the cluster. In their analysis of Abell 2029, they constructed histograms from the 
pixels in radial bins and demonstrated that for each bin the position of the Gaussian peak 
in the pixel intensity histogram is a fair indicator of the ICL levels. The skewness of the 
Gaussian, mainly caused by diffuse light from wings of astronomical objects, is a potential 
source of bias in this measurement. Uson et al. (1991) report that the bias introduced by 
this contamination is about 0.2% of the sky level at one core radius of the cluster if objects 
are left unmasked. As discussed in §2.41 our "5 x" masking scheme reduces the effect below 
the 0.02% level of the sky 

To illustrate our measurement scheme, we show in Figure [5] the histograms that we 
obtain in F775W after masking objects and subtracting stars. We set up radial bins out to 
r ~ 100" (~ 530 kpc) at an interval of Ar = 8". The bottom curve represents the data 
in the r = 8" bin. We displace the subsequent histograms vertically to improve readability. 
The top curve corresponds to the data in the r = 104" bin. The filled squares and thick solid 
line trace the location of the Gaussian peaks. The uncertainties in determining the centroid 
of the Gaussian curve after including photon noise and fiatfielding errors are ~ 0.5 counts 
or /i ~ 30 mag arcsec"^, smaller than the size of the filled square symbol. We choose the 
geometric center of the dark matter ring as the origin of the radial bin. The geometric center 



^The data are available at littp:xnim.esac.esa.int 
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of the dark matter ring is ~ 10"offset toward south from the peak of the X-ray emission (also 
the location of one of the brightest cluster galaxies). 

Figure M shows ten ICL profiles that we measure from both ACS and Subaru data for 
all available passband data. The thick solid lines represent the surface brightness while 
the dashed (dotted) lines show the l-cr limits without (with) including background level 
measurement errors. In case of ACS, we arbitrarily assume that the error in background 
level measurement is at the 0.1% of the sky level (the statistical error is negligible). The 
Subaru profiles have much larger errors mainly because the sky levels are higher. Especially, 
in the z' image the sky level is nearly 3 magnitudes brighter than in the ACS F850LP image, 
and thus the 0.05% flatfielding error is translated into /i ~ 27.5 mag arcsec"^ whereas the 
0.1% flatfielding error in the ACS F850LP corresponds to /i ~ 30 mag arcsec"^. 

Despite the filter-by-filter variation in significance, it is remarkable that all profiles in 
Figure [6] share some conspicuous common features. The surface brightness in every profile 
decreases rapidly from the cluster center to r = 40" — 50". Then, it maintains its level out to 
r ~ 70", where it starts to decrease again for increasing radius. The changes of the slopes at 
the two locations are nearly discontinuous, somewhat more sudden than the corresponding 
changes in mass slopes of J07; the significance of the feature should be evaluated using 
the inner error lines (dashed), which exclude the background level measurement errors. 
This striking consistency across different filters and instruments is the evidence against the 
possibility that the features come from any residual systematic errors. 

The comparison of the ACS and Subaru results at the large radii help us to assess the 
degree of the ICL contamination in the background level for ACS. Although the difference 
in throughput curve between the instruments complicates the detailed comparison, there is 
no apparent indication that any significant systematic underestimation of the ICL level (or 
overestimation of the background level) is introduced in the ACS results; the ICL levels in 
corresponding filters between the two instruments are consistent within the errors. 



3.2. Comparison of ICL Profile to Mass Profile 

The peculiar mass density profile of CI 0024+17 in J07 is characterized by the steep 
decline out to r ~ 50" (~ 265 kpc), the slow rise from r ~ 50" (~ 265 kpc) to ~ 75" 
(~ 397 kpc), and the modest decrease at r > 75". The turnaround at r ~ 75" (~ 397 
kpc) appears as the dark matter ring in the two-dimensional mass reconstruction. The 
ICL profiles in Figure [6] clearly signals these features of the J07 mass profile. Although 
the significance varies across filters and instruments, it is evident that all profiles change 
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their slopes conspicuously at both r ~ 50" and r ~ 72" in the similar way as in the mass 
profile. Considering the size of the uncertainty (we exclude the background level uncertainty 
here because it does not affect the significance of the relative variation), the change of the 
profile slopes is significant in every panel except for the ACS F435W data, for which the 
major source uncertainy comes from a low ICL-to-sky ratio. Of course, the most remarkable 
aspect of the results is the observation that all profiles closely resemble one another, which 
strongly rules out the possibility that the source of these features originate from uncorrected 
instrumental systematics. 

In Figure[7]we overplot these ICL profiles on top of the mass profile of J07 after rescaling 
of the flux via the following form: I'{r) = axl[r)+b. If intracluster stars are truly poissonian 
tracers of underlying dark matter, a strict comparison would require no intercept in the 
transformation. However, considering the uncertainty of the background level determination 
in ICL measurement, the uncertainty of the mass profile caused by mass-sheet degeneracy, 
and the physical mechanism that ICL might be more centrally concentrated, we allow the 
intercept to vary in the fitting procedure. Furthermore, we limit the fitting range to the 
20" < r < 75" region because at r < 20" the ICL measurement is hampered by severe 
masking of bright cluster galaxies, and at r > 75" the ACS data become progressively 
incomplete toward large radii. Also, for some filters ICL profile measurements are sensitive 
to the background estimation in this regime. The black solid line is the mass profile of 
J07. The open circles denote the results from the different filters in ACS whereas the open 
diamonds represent the case for the Subaru images. The thick purple line delineates the 
mean ICL profile. 

This direct comparison more graphically supports the above claim that the ICL of 
CI 0024+17 mimics the peculiar mass profile. The mean ICL profile shows a sudden change 
of slope at r ~ 50" and r ~ 72" as does the mass profile. In addition, the difference among 
the results from different filters is remarkably small. Because the agreement is excellent over 
a large range of radii, this consistency should not be attributed to the rescaling of the profile. 

We note however that the agreement between the ICL and mass profiles degrades at 
r > 72". Apparently, the ICL profile changes its slope more suddenly at r ~ 72" than 
the mass. The difference among the ICL profiles is also slightly larger in this regime. It 
is tempting to view this as indicating that the ICL distribution is more centrally peaked 
than the cluster mass distribution. However, one should use a caution in interpreting the 
difference. The ACS data do not provide complete azimuthal statistics at r ~ 80". Also, the 
ICL level in this range is sensitive to the background level determination. 

The dashed line is an NFW fit to the mass profile using the value in the 20" < r < 50" 
range. Both the cluster mass and ICL profiles deviate from this analytic description, which 
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is frequently used to model halos of relaxed clusters. 



3.3. Is the Peculiar ICL Structure Local? 

The consistent detection of the peculiar ICL feature of CI 0024+17 from both ACS and 
Subaru data provides evidence that the feature does not come from any residual calibration 
errors. And the similarity between the ICL and mass profiles adds further to the credibility 
of the detection. However, one can raise a concern that our masking might miss a few 
low surface brightness objects that fortuitously at a similar distance from the cluster center. 
While we consider that this is very unlikely given the significance and the scale of the features, 
it is still useful to devise a test demonstrating that the features indeed come from a two- 
dimensional, quasi-symmetric diffuse light distribution, not from any localized concentration. 

The most straightforward proof would be the revelation of the discussed feature in the 
two-dimensional ICL map. However, given the low contrast of the structure, it is extremely 
difficult, if not impossible, to produce a convincing diffuse light map of the cluster. The pixel 
histograms in Figure [5] that we discussed in §3.11 help us to quantitatively understand why it 
is next to impossible to see the feature in a two-dimensional map. As discussed already, the 
location of the Gaussian curve centroid is the indicator of the ICL level. The measurement 
leads to an uncertainty small enough (~ 30 mag arcsec"^) for us to give high significance 
to the relative variation of the centroids as a function of radius. However, the width of the 
Gaussian curve itself (i.e., the intrinsic pixel intensity distribution) is much larger than the 
relative variation of the centroids (i.e., ~ 70 counts versus ~ 2 counts), which makes it futile 
to attempt to detect the feature in a two-dimensional map. 

A visual inspection of the two-dimensional ICL map (the lower right panel of Figured]), 
nevertheless, reveals that there is some localized diffuse light in the northwestern quadrant; it 
appears that the diffuse light in this region is associated with several galaxies approximately 
on the r ~ 72" circle. Hence, it is important to examine how significantly the ICL profile 
that we presented in §3.11 is affected by this apparent substructure. To address this issue, 
we divided the cluster field into four azimuthal regions and measured an ICL profile in each 
quadrant. The result from the F775W image is presented in Figure El Although there 
exists a quadrant-to-quadrant variation, the features seen in Figure |6] is also observed in 
each quadrant. Most of all, the ICL level stops decreasing at r = 40" ~ 50" and increases 
out to r ~ 70" in each panel. The sudden change of slope at r ~ 72" is the strongest in 
the NW and SE quadrants. The profile from the NE quadrant shows a weak change of 
the slope at this location whereas the slow rise of the profile between 40 < r < 70" is still 
clear. The SW ICL profile possesses a sudden change of slope at r ~ 50" and the surface 
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brightness is increasing out to r ~ 65". For this quadrant, the change of slope at r ~ 72" 
is subtle, which suggests that the ringlike structure might not be strong there. Nonetheless, 
the interpretation is difficult because the surface brightness level at this location is very low 
{fi > 29 mag arcsec"^) and thus the error is dominated by systematical. 

Therefore, from the above experiment we conclude that the ICL proffie structures re- 
sembling the peculiar mass profile in J07 do not come from any localized substructure, but 
from a two-dimensional, quasi-symmetric distribution of diffuse light in CI 0024-1-17. 



3.4. ICL Color versus Galaxy Color 

The color of ICL relative to the color of cluster galaxies potentially constrains the epoch 
at which most of the intracluster stars were stripped from the cluster galaxies. The literature 
roughly converges on the observation that ICL colors are approximately consistent with the 
cluster galaxy colors, which suggests that intracluster stars may be just as old as cluster 
galaxies. 

Considering the location of the 4000 A break at z = 0.4 in the observed frame, we find 
that a F475W-F625W color gives good contrast between the cluster red sequence and the 
rest of the population in the field. The color magnitude diagram in Figure [H^ with this ffiter 
combination clearly shows the red sequence of the cluster in the ACS field (the red solid line 
approximately delineates the location). We display in Figure [9]d the F475W-F625W color 
profile of the ICL. At large radii, the accuracy highly depends on how well we can determine 
the net amount of ICL. Therefore, the biases in the background level determination cannot 
be neglected here. We assumed that this bias is at the 0.1% of the sky level. 

We find that the ICL color at small radii (r < 40) is consistent with that of the cluster 
red sequence while the color profile tends to become bluer with increasing radius. This trend 
is in agreement with the general expectation that cluster galaxies are bluer at larger cluster- 
centric radii. Robust interpretation of Figure [Hb would be possible if we use a population 
synthesis model assuming several different star formation histories. While we acknowledge 
the importance of this study in determining the age of the progenitor population of the 
ICL, this extensive analysis is beyond the scope of the current paper. However, one possible 
interpretation of Figure [H] is that for the CI 0024-1-17 cluster the observed ICL might be the 
result of on-going stripping of stars (e.g., Moore et al. 1999). If dominant fraction of the 



^The total amount of the ICL flux becomes negative when we subtract the background level. Therefore, 
we arbitrarily adjusted the background value to prevent it 
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ICL stars were stripped a few Gyrs earlier than the current epoch and have evolved passively 
since then, we would expect to observe that the intracluster stars are redder than the cluster 
galaxies because the cluster galaxies continue to form new stars (e.g., Sommer-Larsen et al. 
2005). 



3.5. ICL Fraction 

Together with ICL color, the total amount of light in intracluster stars relative to the 
total cluster light is also an important quantity in helping us to infer the production history of 
intracluster stars. In order to identify the total cluster light, we combined the photometric 
redshift catalog of J07, and the publicly available spectroscopic catalog of Moran et al. 
(2005). Figure 3 and 4 of J07 demonstrate that the photometric redshifts estimated from six 
ACS filters with Hubble Deep Field North (HDFN) prior give highly consistent results with 
the spectroscopic measurements and reliably identify the cluster members. We replaced the 
photometric redshifts if their spectroscopic redshifts are known (~ 140 objects) and used 
this merged catalog to select the cluster members of CI 0024+17. The rest of the procedure 
is identical to the one in our ICL measurement described in §3. II except that 1) we did not 
mask out the cluster members and 2) took the mean value in each annulus instead of the 
location of the Gaussian peak. 

Figure [10] shows the ICL fraction measured from the five ACS filter images; we omit the 
results for F435W, whose error bars are substantially larger than the rest. As in the case for 
the ICL color measurement, the ICL fraction at large radii is sensitive to the background level 
determination and we again assume 0.1% of the sky level for the uncertainty. It appears 
that the mean level of the ICL fraction somewhat depends on the used filter. Although 
this observed trend possesses low significance at large radii due to the systematics in the 
background level determination, the results at small radii (r < 40") show that the ICL 
fraction is on average higher in blue filters (i..e., F475W, F555W, and F625W). 

In order to estimate the cumulative ICL fraction within an aperture, we need to assume 
the behavior of the ICL profile at r < 15", where our masking procedure left no area for direct 
ICL measurement. Extrapolating the central value using the trend at r < 40", we obtain the 
results in Table 4. The uncertainties are determined from Monte-Carlo simulations while the 
correlations between points are taken into account. The literature does not converge on the 
ICL fraction, and perhaps this represents rather a wide range of values for different clusters 
than scatters in measurements. Certainly, the ICL fraction (> 30%) of CI 0024+17 that we 
quote here indicates that the cluster might be near the high end of the distribution. Tyson et 
al. (1998) determined the ICL fraction for CI 0024+17 to be 15%±3% within the rsimlOO kpc 
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region. Because their analysis is based on the Wide Field Planetary Camera 2 (WFPC2) 
data, which covers only ~ 50 % of the ACS field, we suspect that their background level 
measured within the WFPC2 field might be biased high due to the significant ICL level; the 
ICL profile in our study shows that the surface brightness level remains high out to r ~ 70". 

Given the numerical studies suggesting that most of intracluster stars become unbound 
at ^ < 1 with no preferred epoch (e.g., Murante et al. 2007), the high ICL fraction of 
CI 0024+17 at z = 0.4 is somewhat unexpected. However, if the high-speed encounter 
scenario of Czoskey et al. (2001) for CI 0024+17 is considered, it is plausible that the 
collision of two equal-mass subclusters 1-2 Gyrs ago might have liberated unusual amount 
of stars in the cluster galaxies. 



4. INTRACLUSTER MEDIUM PROFILE OF CI 0024+17 

Although the coUisional and dissipative nature makes its spatial distribution often differ 
from that of dark matter, intracluster gas distribution in a galaxy cluster strongly signals 
the structure of the cluster potential dominantly determined by the cluster dark matter. 
Therefore, whether or not one can find a good spatial agreement between the two results, a 
detailed comparison often provides an insight into the dynamical history of the two cluster 
constituents (e.g., Jee et al. 2005; Clowe et al. 2007). 

Our previous analysis of the Chandra data of CI 0024+17 confirmed the unusual X-ray 
profile of the cluster originally reported by Ota et al. (2005). The central excess of the X-ray 
photons only allows us to describe the profile with a superposition of two isothermal profiles 
of different (3 parameters. Nevertheless, the X-ray surface brightness profile of J07 did not 
show any features reminiscent of the reported dark matter profile. 

However, repeating the J07 analysis with a different choice of the cluster center shows 
that there presents a 3.7cr bump in the X-ray profile at r ~ 60" as shown in Figure [TTJ The 
result in the left panel is obtained when the center of the dark matter ring is chosen as the 
origin of the radial bins whereas the the right panel shows the X-ray profile when the center 
is placed on the peak of the X-ray emission as in the J07 analysis. It is worth noting that 
the distance between these two locations is only ~ 10". The fact that this small change 
greatly alters the significance of the feature indicates that the contrast is low and the bump 
is associated not with any unidentified point source, but with smooth, azimuthal distribution 
of excess X-ray photons. 

The X-ray surface brightness stops decreasing at r ~ 50" and increases slowly until it 
suddenly drops again at r ~ 60". The overall behavior is reminiscent of the features seen in 
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both the ICL and mass profiles of the cluster. 

Similar to the case in our ICL measurement, this X-ray feature becomes more credible if 
a consistent structure is observed in data from a different instrument. We present the result 
from the XMM-Newton in Figure [TTJ The large PSF of the instrument certainly smooths 
the aforementioned feature in the X-ray profile. However, the profile still shows a 2.4 a bump 
at the same location. 



5. DISCUSSION 

5.1. Discrepant Mass Profiles in the Literature 

The galaxy cluster CI 0024+17 has been known for its large (a factor of two or more) 
mass discrepancy between gravitational lensing and X-ray results. However, the large mass 
discrepancy from different lensing analyses has not received its due attention. We compare 
five different mass profiles found in the literature in Figure [131 Although some of the large 
discrepancy at small radii (< 20") is due to the difference in the choice of the cluster center, 
it is remarkable that different profiles give a very wide range of mass density outside the 
location of the critical curve (~ 30"). 

Some of the discrepancy might be explained by mass-sheet degeneracy in gravitational 
lensing. Mass-sheet degeneracy refers to the invariant lensing observable under the following 
transformation of the mass density: k ^ Xn + {1 — A). When k, is small (<< 1), the 
transformation can be viewed as adding an additional sheet of mass (hence, the term "mass- 
sheet "). When K is not small, the transformation becomes more conspicuous in change of 
slope. Indeed, most of the different mass profiles in Figure [13] roughly overlap under this 
transformation with a proper choice of A and a centroid. 

However, the J07 profile is unique in the sense that the density does not decrease at 
r > 50". As repeatedly discussed above, the mass density increases from this point on until 
it drops at r ~ 75". Therefore, it would be incorrect to claim that the peculiarity of the 
J07 mass profile is attributed to unresolved mass-sheet degeneracy because the invariant 
transformation does not change the sign of the slope! 

If the multiple image identification by Zitrin et al. (2009) is flawless and the bias in 
their photometric redshift estimation is negligible, their mass profile should be regarded as 
the most accurately calibrated profile to date because the 33 multiply-imaged sources at 
different redshifts certainly provide strong constraints on the mass slope. Nevertheless, it is 
important to remember that these strong-lensing features are not useful in constraining the 
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mass profile at r > 50", where no multiple images have been identified. Hence, it may be 
possible that the J07 mass profile would be transformed to match the Zitrin et al. (2009) 
profile at r < 50". Nevertheless, the result should still reveal the peculiar structure at r > 50" 
(with reduced density because A must be greater than unity in this case). A further analysis 
combining the strong lensing constraints of Zitrin et al. (2009) with the weak-lensing data 
is a subject of our future investigation (M. Jee et al. in prep). 

The diverse results in the literature shown in Figure [12] clearly demonstrates that the 
mass structure of the cluster CI 0024+17 is still disputed in the community, and thus more 
concerted efforts toward a better mass model, such as spectroscopic redshift survey of the 
multiply-lensed sources, high-resolution infrared imaging aimed for better photometric red- 
shift estimation of background galaxies, etc., are desired. On the other hand, it is important 
to investigate the issue with different observables. Our study of the intracluster stars as 
potential tracers of the underlying dark matter structure is one such approach, and here we 
argue that the peculiar features of the ICL that we report in this paper strongly favors the 
J07 result among the five models. Of course, it is still an open question how well the intra- 
cluster stars trace the underlying dark matter. While there is a consensus that intracluster 
stars are not bound to individual galaxies but to the cluster potential, the answer to this 
question certainly depends on when, where, and how the ICL is produced. We discuss the 
issue in ^5.2[ 



5.2. Interpretation of the Similarity between ICL and Mass Profiles 

If intracluster stars are dominantly produced by truncation of galaxy halos during the 
initial formation of a cluster (e.g., Merritt 1984), the dynamical history of the intraclus- 
ter stars is as old as the cluster itself. Hence, in this scenario we expect to observe the 
strongest correlation between ICL and mass distributions. Many A^-body simulations, how- 
ever, indicate that other mechanisms such as high-speed encounters between cluster galaxies 
(Moore et al. 1996), stripping from a cluster potential (Byrd & Valtonen 1990), stripping in 
galaxy groups falling into the cluster (Mihos 2004), etc., might also contribute significantly 
to the overall production, which then would make it less straightforward to interpret the ICL 
structure in conjunction with the mass distribution. It is an open question which of these 
mechanisms is most dominant. 

Recently, Murante et al. (2007) investigated the origin of the intracluster stars through a 
cosmological hydrodynamical simulation. They reported that a dominant fraction (~ 50%) of 
the intracluster stars are produced in the formation of the most massive galaxy at the cluster 
center, which is consistent with the observation that the ICL is more centrally concentrated 
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than the cluster galaxies (e.g., Zibetti et al. 2005). If cluster galaxies are poissonian, albeit 
sparse, tracers of underlying dark matter, this suggests that perhaps the ICL profile might 
be more peaked also than the cluster mass profile. 

With this view in mind, the resemblance of the ICL profile to the mass profile shown 
in Figure [7] may seem somewhat surprising. However, it is important to remember that 
the ICL profiles are rescaled to match the mass profile, and thus it is trivially obvious that 
the agreement in the range 20" ^ r < 50" indicates nothing more than that both the ICL 
and mass profiles monotonously decreases in this regime. What deserves our attention, 
nonetheless, is that both profiles change the slope at r ~ 50" by the same amount (of course, 
after rescaling is applied), giving a good agreement out to r ~ 70". The average ICL levels 
between the r ~ 20" and r ~ 50" regions differ by a factor of five while in mass the density 
ratio is less than a factor of two. Because the mean slope at 50 < r < 70" is much smaller 
than the value at r < 50", the ICL profile, although being much steeper than the mass 
profile, could be brought to agree after rescaling. 

Therefore, our results presented here should not be interpreted as a demonstration that 
the intracluster stars in CI 0024+17 can be used as direct (i.e., poissonian) tracers of the 
underlying dark matter. Our ICL profile measurement supports the view that the dominant 
ICL production mechanism is associated with the formation of the brightest cluster galaxies 
at the cluster center, and thus their distribution is expected to be more centrally concentrated 
on average. However, dynamically old intracluster stars are given the chance to travel further 
away from the cluster center and mix with the cluster halo dark matter. Hence, the ICL 
profile should represent the dynamical history of the cluster in the central region, as well as 
the ICL production history. 

The expectation that cluster outskirts might be severely depleted of intracluster stars 
may aid us to explain the difference between the mass and ICL profiles at r > 72". As 
already discussed in §3.2[ the ICL profiles change their slopes more suddenly at r ~ 72" 
than the mass profile, and hereupon the ICL levels quickly approach the background level. 
It is plausible that this is because the initial ICL level is already low in this region before the 
dark matter structure is disrupted. Alternatively, we can consider the possibility that the 
J07 mass profile is biased at r > 75", and the real mass profile is similar to the level that is 
indicated by the ICL profiles. The radial mass profile in J07 does not represent the complete 
azimuthal average at r > 85", and thus the employed regularization can progressively bias 
the mass density high as r increases. 
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5.3. Is the Cluster Galaxy Distribution Inconsistent with the Mass Profile of 

JOT? 

The current, well-accepted, structure formation paradigm postulates that galaxies are 
biased tracers of hosting dark matter. Gravitational lensing studies of galaxy clusters have 
demonstrated that even in clusters of complex morphology the spatial correlation between 
galaxies and dark matter is strong (e.g., Hoekstra et al. 2000; Jee et al. 2005; Clowe et al. 
2006), although Mahdavi et al. (2007) recently claims that a puzzling counter example is 
seen in Abell 520. 

If galaxies and dark matter indeed follow each other, can we detect the ringlike structure 
of J07 in the cluster galaxy distribution? Of course, given that the ringlike structure of J07 
has an only 5% contrast, it is impossible to detect the feature as an enhancement of galaxy 
number density around the "ring" even if we improve the statistics by azimuthally averaging 
the numbers. Figure [H] shows the projected number density of the cluster members as a 
function of radius (0.37 < z < 0.41) obtained from the publicly available catalog of Moran 
et al. (2005), which extends the Czoske et al. (2001) catalog. The shot noise at r ~ 75" 
is ~ 45%, which is about 9 times higher than the contrast; we did not correct for the 
completeness of the spectroscopic survey, which is negligible in this region compared with 
the shot noise. 

Apart from the issue of detecting the ringlike structure of J07, it is still interesting 
to examine if the overall galaxy number density profile is consistent with the mass profile. 
As mentioned by J07, the mass profile at r > 80" is measured from an incomplete circle. 
Therefore, in order to enable a fair comparison, we should limit the fitting range to r < 80". 
The solid line in Figure [T3] shows that the J07 mass profile is consistent with the galaxy 
distribution {x^/d.o.f. = 0.45) within r < 80". When we instead try the fit with the data 
points at r < 100", the result (dotted) is severely affected by the data at r > 80", giving 
X^/d.o.f. = 1.92. Qin et al. (2008) used the latter comparison to argue that the galaxy 
distribution is inconsistent with the J07 mass profile. Because the mass profile of J07 at 
r > 80" does not represent the azimuthal average, we do not agree with their interpretation. 

How many galaxies do we need to detect the peculiar dark matter structure? Within 
the r = 100" (80") circle, there are 99 (84) 0.37 < z < 0.42 galaxies in the Moran et al. 
(2005) catalog. The mass density at the r ~ 75" bump is about 5% higher than the value at 
r ~ 50". If we want to detect this rising mass density in the 50" < r < 75" region with the 
galaxy distribution at the 2 a level, we need ~ 50 times more galaxies, which is even greater 
than the total number of objects (~ 4000) detected in the ACS image! Our photo-z analysis 
estimates that about 350 out of 4000 galaxies are in the redshift range 0.3 < z < 0.5, and 
this number certainly does not help us to reveal the ringlike dark matter structure. 
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6. CONCLUSIONS 

Our main results of the investigation of the ICL and ICM in the galaxy cluster CI 0024+17 
are summarized as follows. 

• The ICL profile closely resembles the peculiar dark matter structure reported in J07, 
which stops decreasing at r ~ 50" (~ 265 kpc) and slowly increases until it turns over 
at r-75" (- 397 kpc). 

• The feature is present in both ACS and Subaru images in nearly all available passband 
images with a significance in general higher in red filters. 

• The ICM profile shows a bump at r ~ 60" (~ 318 kpc). The radius of the ICM ring is 
~ 20% smaller than that of the dark matter ring. 



• 



• 



• 



The two dimensional map of the diffuse light in the cluster shows that the ICL distri- 
bution is elongated toward Northwest following cluster members. 

The ICL color is consistent with the cluster red-sequence color at small radii (r < 40"), 
but becomes bluer for increasing radii. 

The ICL fraction of CI 0024+17 is ~ 28% (~ 37%) in F775W and F850LP (F475W, 
F555W, F625W) without any strong indication of radial dependence. Considering the 
redshift of the cluster, these values are likely to be at the high end of the distribution. 

These results not only support the presence of the peculiar dark matter structure of 
J07 revealed by coherent fluctuation of the background galaxy shapes across the r ~ 75" 
circle, but also demonstrates the usefulness of the ICL as a powerful probe of dark matter 
substructure in galaxy clusters. 

Theoretically, however, it is still disputable whether or not the ICL distribution for a 
z = 0.4 cluster should reveal the underlying dark matter. Numerical simulations by Rudick 
et al. (2006) show that the ICL morphology is a strong function of time. At early epoch most 
intracluster stars reside in outer halos of individual galaxies while with time more tidal and 
fllamentary structures develop. As the evolution progresses further, it is shown that these 
transitory features disappear into the more diffuse, common cluster halo. If we associate the 
redshift z = OA with the young dynamical age in ICL evolution, it is difficult to understand 
why the cluster's ICL follows the dark matter structure so remarkably well as shown in this 
work. On the other hand, if CI 0024+17 is regarded dynamically already mature based on its 
optical richness or mass, the similarity between the ICL and dark matter structures may not 
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be surprising. The high ICL fraction of the cluster supports the possibihty that CI 0024+17 
is aheady mature in the context of the ICL evolution. Alternatively, it is also plausible that 
the high-speed collision 1-2 Gyrs ago might have liberated unusual amount of stars from the 
cluster galaxies. 

The detection of the ICM bump at r ~ 60" opens an important opportunity to study 
the dynamical interplay between the intracluster medium and the underlying dark matter. 
One simplistic scenario that one can make out of this feature is that both the ICM and dark 
matter ringlike structures have the same dynamical origin, and both rings (or shells) are 
expanding. In this case, the expansion of the ICM ring must lag behind the dark matter 
ring because of the ram pressure, which may explain the observed ~ 15" (~ 80 kpc) offset. 

In a recent dark matter only simulation, ZuHone et al. (2008) claim that a bump can 
arise in cluster collisions only when the initial particle velocity distribution is circular, and 
the feature becomes a "shoulder" as the tangential anisotropy decreases. Our re-examination 
of the J07 collision simulation reveals that the initial orbits of the particles in the cluster 
core were indeed biased toward tangential anisotropy although the velocity anisotropy is not 
as extreme (i.e., circular) as shown by ZuHone et al. (2008). As the distinction between 
"bump" and "shoulder" is whether or not the slope of the post-collision mass profile in the 
intermediate range (i.e., 50" < r < 75") is positive, we believe that many other parameters 
such as the slope of the initial density profile, the fraction of the gas particles, the relative 
mass ratio/speed of the two subclusters, etc. should be also considered together with the 
initial distribution of the particle orbits in order to fully address the issue. Nevertheless, it 
is clear from the ZuHone et al. (2008) simulation that the creation of the ringlike structure, 
if the high-speed collision is to be the cause as suggested by J07, requires some degree of 
tangential anisotropy at least in the core of the cluster. 

We acknowledge Maxim Markevitch for first informing us of the presence of the ringlike 
structure in the Chandra image of CI 0024+17. We would like to thank Bernard Fort for 
suggesting this study at the 23rd lAP colloquium. Furthermore, we would like to thank 
Tony Tyson, Holland Ford, David Wittman, and John Blakeslee for useful discussions. This 
research is in part supported by the grant from the TABASCO foundation awarded in the 
form of the Large Synoptic Survey Telescope Cosmology Fellowship. 
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ACS F475W+F625W+F850LP ■ Diffuse Light in A 



Fig. 1. — ACS and Subaru color images of CI 0024+17 and two-dimensional ICL map. The 
top panel shows the central 10' x 10' region of the color-composite Subaru image. White 
circles denote spectroscopic cluster members (0.37 < z < 0.42). The yellow outline illustrates 
the 3' X 3' region observed with ACS. The ACS color-composite image is shown separately 
in the lower left panel. The two red dashed-circles delineate the r ~ 40" and r ~ 70" radii, 
across which a sudden change of slope in the ICL profile occurs. The lower right panel 
displays the diffuse light in ACS F625W image. This two-dimensional ICL image is created 
by masking out astronomical objects, subtracting bright stars, and then applying 3" x 3" 
box median-smoothing. The color scale is linear. The diffuse light extends towards the 
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0.97 



1.02 




Fig. 2. — Residual sky flats in the six ACS/WFC broadband filters. These flats are created 
from the pipeline-reduced FLT files, and should reveal any residual features that are not 
corrected by the pipeline LP-flats (see text for details in the procedure). The displayed gray 
scale is linear. The rms residual is ~ 0.5%. The most outstanding large scale feature is the 
donut-like pattern particularly clear in F775W and F850LP. This residual feature may arise 
from the limited accuracy of the polynomial interpolation in the L-fiat creation or from the 
difference of the spectra between the sky and the 47 Tuc stars as suggested by Pavlovsky 
et al. (2005). One of the most conspicuous small scale features is the gridlike pattern. We 
verified that this is not an artifact of the 32x32 box median smoothing that we apply before 
stacking is carried out. Although we suspect that the feature might be associated with the 
CCD fabrication procedure, the exact cause of this ~ 60 pixel-scale feature is still under 
investigation. In addition, the residual flats reveal some new dust-moats and under- sensitive 
CCD columns. 
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Fig. 3. — Experiment with different masking sizes. We investigated the effect of the masking 
size on the centroids of the Gaussian curves by varying masking sizes and repeating the 
measurements. As an example, shown here is the result from the r ~ 80" bin in F775W. The 
different Gaussian curves represent the results from different masking schemes. The vertical 
dashed lines show the corresponding centroids. We observe that the centroid converges when 
the major and minor axes of the masking ellipse are three times or greater than SExtractor 
measurements. We conservatively selected five times the values given by SExtractor in 
defining our masking ellipses throughout the analysis. 
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Fig. 4. — Extended PSF wings of ACS/WFC (left) and Subaru/Suprime-cam (right). The 
normahzation is chosen in an arbitrary way so that the profiles from different filters approx- 
imately overlap at r < 20". The ACS results are obtained from the archival observation of 
the star HD39060 (Beta Pictoris). The Subaru PSFs were measured from the brightest star 
(r ~ 8.6 mag) in the CI 0024+17 image (~ 7') away from the cluster center. 
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Fig. 5. — Pixel intensity distribution in the ACS/WFC F775W image. Objects were masked 
and bright stars were subtracted before the statistics was evaluated. Each curve represents 
the pixel intensity distribution in each radial bin. The bottom curve represents the data in 
the r = 8" bin. To improve readability, we displace the histograms vertically. The top curve 
corresponds to the data in the r = 104" bin. We adopt the center of the Gaussian peak as 
the ICL level of the bin (filled square and thick solid line). The uncertainties in determining 
the centroid of the Gaussian curves after including photon noise and flatfielding errors are 
~ 0.5 counts, less than the size of the filled square symbols. 
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Fig. 6. — ICL profile in CI 0024+17. The thick sohd hnes represent the surface brightness 
while the dashed (dotted) lines show the 1-a limits without (with) errors in background level 
determination included. In case of ACS, we arbitrarily assume that the error in background 
level measurement is at the 0.1% of the sky level. Because of the limited field of view, 
ICL measurement from the ACS data is obtained from an incomplete circle at r > 85" 
(dot-dashed). 
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Fig. 7. — Comparison of the ICL profiles with the mass density profile. The ICL profiles 
from nearly all the filters reveal the critical features of the mass profile: the steep decline at 
r < 50", the slow rise at 50" ^ t < 75", and the turnaround at r ~ 75". However, we note 
that the ICL profiles decrease faster after their turnover at r ~ 72" than the mass profile 
of J07. We display the NFW fit result (obtained by forcing the profile to match the mass 
at 20" < r < 40") in order to illustrate how much both the cluster mass and ICL profiles 
deviate from this analytic model frequently used to describe relaxed clusters. 
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Fig. 8. — Same as in Figure |6] except here we measure the ICL profile in each quadrant 
of the ACS F775W image. Despite the quadrant-to-quadrant variations, the signal of the 
peculiar ICL structure is present in all four quadrants. The dashed-line does not include the 
uncertainty from background measurement. 
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Fig. 9. — Galaxy color versus ICL color, (a) Color magnitude diagram in the ACS field. 
We used SExtractor's isophotal magnitude to compute the galaxy colors while the F625W 
magnitude in the X-axis is SExtractor's MAG_AUTO. The red solid line approximately 
indicates the location of the cluster red-sequence, (b) ICL color profile as a function of 
radius. The dotted lines represent the 1-a uncertainty including errors in photon statistics, 
flatfielding inaccuracy, and background level determination. We assume that the error in 
the background level determination is at the 0.1% level of the sky. 
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Fig. 10. — Ratio of ICL flux to total cluster flux. We combined the photometric redshift 
catalog of JOT, and the publicly available spectroscopic catalog of Moran et al. (2005) to 
identify the cluster members. 
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Fig. 11. — X-ray surface brightness profile created from the Chandra data, (a) When the 
origin agrees with the center of the dark matter ring, it reveals an overdensity at r ~ 60" 
(the inset picture shows the point-source removed X-ray image and the '+' symbol denotes 
the origin of the radial bin), (b) When we align the origin to the peak of the X-ray emission 
(also the location of the mass and luminosity peak), the bump in (a) becomes much less 
prominent. We let the neighboring bins overlap by ~ 50% to suppress the shot noise. 
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Fig. 12. — X-ray surface brightness profile from the XMM-Newton data. The bump at 
r ~ 60" seen in Figure [11] is somewhat smoothed by the larger PSF but still visible as a 
shoulder. The feature is difficult to observe if we chose the X-ray peak as the origin of the 
radial bin as in the case of Figure [Til The scale is arbitrary. The neighboring bins do not 
overlap. 
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Fig. 13. — Discrepant mass profiles of CI 0024+17 in tlie literature. The unit is the critical 
lensing mass density at the z = 3 reference redshift. For the Broadhurst et al. (2000) results, 
we reproduce their two-dimensional mass map, and measure the profile by placing the center 
at the geometric center of the ringlike structure of J07. For the Tyson et al. (2000) and 
Kneib et al. (2003) we use their parametric fit to the mass profiles. The Zitrin et al. (2009) 
profile shown here is their best-fitting model (read off from Figure 5 of the paper). 
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Fig. 14. — Projected galaxy number density versus radius. We use the publicly available 
spectroscopic catalog of Moran et al. (2005) to plot the number density of the spectroscopic 
members in the cluster as a function of radius (filled circle). Because the mass profile of 
JOT is progressively biased beyond r > 80" (dashed line) because of the limited field of 
view, a fair comparison may require limiting the range to fit at r < 80". This gives a 
goodness-of-the-fit of x^/d.o.f = 0.45 (solid), which illustrates that the galaxy distribution 
is statistically consistent with the J07 mass profile. If the data points at r > 80" are included, 
the goodness-of-the-fit degrades to x^ /d.o.f = 1.92 (dotted). However, interpretation needs 
a caution because the result is severely affected by the data points at r > 80", which have 
smaller error bars yet with lower completeness. 
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Table 1. ACS Image Sky Levels 



Filter 



Visits 



Total Exp. Time 
(s) 



Mean Sky Level 
(mag arcsec^'^) 



Sky Level Variation Between Visits 
(%) 



F435W 


5 


6435 


23.1 


F475W 


4 


5072 


22.9 


F555W 


4 


5072 


22.6 


F625W 


7 


8971 


22.2 


F775W 


8 


10144 


22.2 


F850LP 


6 


16328 


22.1 



11 

4 
3 
8 
2 
2 



Table 2. Subaru Image Sky Levels 



Filter 


Visits 


Total 


Exp. Time 


Mean Sky 


Level at zenith 


Sky Level Variation Between Visits 








(s) 


(mag 


arcsec^^) 


(%) 


B 


3 




3,600 




22.6 


2 


Re 


11 




5,280 




21.2 


3 


z' 


8 




1,980 




19.2 


5 


NBgj2 


6 




10,800 




20.7 


11 



Table 3. ACS Sky Flat Accuracy 



Filter 



Number of Used Images 



Deviation from Pipeline Flats Accuracy 

(%) (%) 



F435W 


90 


F475W 


138 


F555W 


92 


F625W 


155 


F775W 


154 


F850LP 


312 



0.5 
0.4 
0.4 
0.5 
0.4 
0.7 



0.08 
0.07 
0.06 
0.07 
0.07 
0.08 



Table 4. Subaru Sky Flat Accuracy 



Filter 


Number of Used Images 


Accuracy 

(%) 


B 


49 


0.08 


Re 


138 


0.05 


z' 


155 


0.05 


NBgi2 


41 


0.09 
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Table 5. Errors in Subaru Background Level Estimation 



Filter 


Fraction of 


Sky 


Surface Brightness 




(%) 




(mag arcsec"^) 


B 


0.05 




30.9 


Re 


0.04 




29.7 


1 


0.04 




27.6 


NBgr2 


0.08 




28.4 



Table 6. ICL fraction in CI 0024+17 



Filter 


r < 200 kpc 


r < 500 kpc 




(%) 


(%) 


F475W 


35 ±6 


38 ±12 


F555W 


35 ±4 


37 ±9 


F625W 


35 ±2 


35 ±5 


F775W 


27 ±1 


26 ±3 


F850LP 


29 ±1 


30 ±3 
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